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Abstract 



Klee's Measure Problem (KMP) asks for the volume of the union of n axis-aligned boxes in R"^. 
Omitting logarithmic factors, the best algorithm has runtime 0*{n''^'^) [Overmars,Yap'91]. There are 
faster algorithms known for several special cases: Cube-KMP (where all boxes are cubes), Unitcube- 
KMP (where all boxes are cubes of equal side length), Hypervolume (where all boxes share a vertex), 
and fc-GROUNDED (where the projection onto the first k dimensions is a Hypervolume instance). 

In this paper we bring some order to these special cases by providing reductions among them. In 
addition to the trivial inclusions, we establish Hypervolume as the easiest of these special cases, and 
show that the runtimes of Unitcube-KMP and Cube-KMP are polynomially related. More importantly, 
we show that any algorithm for one of the special cases with runtime T{n, d) implies an algorithm for the 
general case with runtime T(n, 2d), yielding the first non-trivial relation between KMP and its special 
cases. This allows to transfer W[l]-hardness of KMP to all special cases, proving that no n"'^'*' algorithm 
exists for any of the special cases under reasonable complexity theoretic assumptions. Furthermore, 
assuming that there is no tmproved algorithm for the general case of KMP (no algorithm with runtime 
0{n'''^^~^)) this reduction shows that there is no algorithm with runtime 0{n^''^'^^ ^^~'^) for any of the 
special cases. Under the same assumption we show a tight lower bound for a recent algorithm for 2- 
Grounded [Yildiz,SuriT2]. 
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1 Introduction 



Klee's measure problem (KMP) asks for the volume of the union of n axis-aligned boxes in M**, where d 
is considered to be a constant. This is a classic problem with a long history [2, 8, 9, 13, 15, 16]. The fastest 
algorithm has runtime C(n''/^ logn) for d ^ 2, given by Overmars and Yap [15], which was slightly improved 
to 7xrf/220(iog* n) Chan [8]. Thus, for over twenty years there has been no improvement over the runtime 
bound n''/^. As already expressed in [8], one might conjecture that no improved algorithm for KMP exists, 
i.e., no algorithm with runtime 0{n'^^^~^) for some e > 0. 

However, no matching lower bound is known, not even under reasonable complexity theoretic assumptions. 
The best unconditional lower bound is O(nlogn) for any dimension d [9]. Moreover, Chan [8] proved that 
KMP is W[l]-hard by giving a reduction to the A;-Clique problem. Since his reduction has k = Q{d), we can 
transfer runtime lower bounds from fc-Clique to KMP, implying that there is no n°^'^^ algorithm for KMP 
under reasonable complcixity theoretic assumptions (assuming the Exponential Time Hypothesis, see [14]). 
However, this does not determine the correct constant in the exponent. 

Since no further progress was made for KMP for a long time, research turned to the study of special 
cases. Over the years, the following special cases have been investigated. As we have no space for a detailled 
listing of the literature on all special cases, we only mention the asymptotically fastest results. 

• Cube-KMP: Here the given boxes are cubes, not necessarily all with the same side length. This 
case can be solved in time 0(71^**+^^/^) for d ^ 2 [6]. In dimension d = 3 this has been improved to 
0{nlog^ n) by Agarwal [1]. In dimensions d ^2 even the general case can be solved in time 0{nlogn), 
the same bound clearly applies to this special case. As described in [6], there are simple reductions 
showing that the case of cubes is roughly the same as the case of "a-fat boxes" , where all side lengths 
of a box differ by at most a constant factor a. 

• Unitcube-KMP: Here the given boxes are cubes, all of the same side length. Since the combinatorial 

complexity of a union of unit cubes is C'(nL''/^J) [5], there is an immediate algorithm with runtime 
0(„L'i/2j polylogn) [12], improving upon the best general case algorithm in odd dimensions. Further- 
more, there is an elaborate algorithm running in time 0(n^'^/^^"^+T3k polylogn) [7], improving m even 
dimensions. Again, there is a generalization to "a-fat boxes of roughly equal size", and any algorithm 
for Unitcube-KMP can be adapted to an algorithm for this generalization [6]. 

• Hypervolume: Here all boxes have a common vertex. Without loss of generality, we can assume 
that they share the vertex (0, . . . , 0) e R*^ and lie in the positive orthant ]R>o- This special case is of 
particular interest for practice, as it is used as an indicator of the quality of a set of points in the field 
of Evolutionary Multi-Objective Optimization [3, 11, 18, 19]. The only asymptotic improvement known 
(compared to the general case of KMP) is an algorithm with runtime 0(71 logn) for d = 3 [4]. The 
same paper also shows an unconditional lower bound of f2(nlogn) for d > 1. 

• fc-Grounded: Here the projection of the input boxes to the first k dimensions is a Hypervolume 

instance, where ^ fc ^ rf, the other coordinates are arbitrary. This rather novel special case appeared 
in [17], where an algorithm with runtime C(n'^''^^^/^ log n) for d^'i was given for 2-GrOUNDED. 

Note that for none of these special cases W[l]-hardness is known, so there is no larger lower bound than 
O(nlogn), not even under reasonable complexity theoretic assumptions. Also note that there are trivial 
inclusions of some of these special cases: Each special case can be seen as a subset of all instances of the 
general case. As such subsets, the following inclusions hold. 

• Unitcube-KMP C Cube-KMP c KMP. 

• {k + 1)-Grounded C /j-Grounded for all k. 

• rf-GROUNDED = Hypervolume and 0-Gr.ounded = KMP. 

Clearly, for such an inclusion it holds that any algorithm for the superset is also an algorithm for the subset, 
with the same runtime. Moreover, any lower bound for the subset is also a lower bound for the superset. 
This allows to transfer some results listed above to other special cases, e.g., the CuBE-KMP algorithm with 
runtime 0{n^'^~^'^^/^) also applies to Unitcube-KMP. 
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1.1 Our results 



We present several reductions among the above four special cases and the general case of KMP. They provide 
bounds on the runtimes needed for these variants and, thus, yield some order among the special cases. 
Our first reduction relates Hypervolume and Unitcube-KMP. 

Theorem 1. If there is an algorithm for Unitcube-KMP with runtime T'uNiTcuBE-KMp(f^) ci); then there is 
an algorithm for Hypervolume with runtime 

^Hypervolume C'^) d) ^ C'(T'unit(;;ube- 

Note that if Hypervolume were a subset of Unitcube-KMP, then the same statement would hold, with 
the constant hidden by the 0-notation being 1. Hence, this reduction can nearly be seen as an inclusion. 
Also note that this reduction allows to transfer runtime bounds from Unitcube-KMP and CuBE-KMP to 
Hypervolume, e.g., there is an algorithm for Hypervolume with runtime 0(n(''+^)/^) for d > 2. Moreover, 
together with the trivial inclusions this reduction establishes HYPERVOLUME as the easiest of all studied 
special cases. In particular, lower bounds for Hypervolume also holds for all other special cases. 

Corollary 2. For all studied special cases, Hypervolume, Unitcube-KMP, Cube-KMP, and k- 
Grounded (for any ^ k ^ d), we have the unconditional lower bound n{n\ogn) for any d > 1. 

One can find contradicting statements regarding the feasibility of a reduction as in Theorem 1 in the 
literature. On the one hand, existence of such a reduction has been mentioned in [17]. On the other hand, a 
newer paper [10] contains this sentence: "Better bounds have been obtained for the KMP on unit cubes 
but reducing the hypervolume indicator to such problems is not possible in general." In any case, a proof of 
such a statement cannot be found anywhere in the literature. 

Our second reduction substantiates the intuition that the special cases Cube-KMP and Unitcube- 
KMP are very similar, by showing that their runtimes differ by at most a factor of 0{n). Recall that 
Unitcube-KMP C Cube-KMP was one of the trivial inclusions. Clearly, it implies that if there is an 
algorithm for CuBE-KMP with runtime T{n,d), then there is an algorithm for Unitcube-KMP with the 
same runtime. We prove an inequality in the other direction in the following theorem. 

Theorem 3. If there is an algorithm for Unitcube-KMP with ranfime Tunitcube-kmp('^, rf); then there is 
an algorithm for Cube-KMP with runtime 

T'cuBE-KMp(?^) d) ^ 0(n ■ Tunitcube-KMP 

Our third and last reduction finally allows to show lower bounds for all special cases. We show an 
inequality between the general case of KMP and 2fc-GROUNDED, in the opposite direction than the trivial 
inclusions. For this, we have to increase the dimension in which we consider 2fc-GROUNDED. 

Theorem 4. // there is an algorithm for 2fc-GROUNDED in dimension d + k with runtime T2fc_GRouNDED d+ 
k), then there is an algorithm for KMP in dimension d with runtime 

TKMp(n, d) ^ 0(T'2fe_GRouNDED(", d + k)). 

Note that, if we set k = d, the special case 2A:-Grounded in d + k dimensions becomes Hypervolume 
in 2d dimensions. Since we established Hypervolume as the easiest variant, the above reduction allows to 

transfer W[l]-hardness from the general case to all special cases. Since the dimension is increased only by a 
constant factor, even the tight lower bound on the runtime can be transfered to all special cases. 

Corollary 5. There is no n"^"^^ algorithm for any of the special cases Hypervolume, Unitcube-KMP, 
Cube-KMP, and A;-Grounded, under the Exponential Time Hypothesis. 

We get more precise lower bounds if we assume that no improved algorithm exists for KMP (no algorithm 
with runtime 0{n'^^^~'^)). Since any algorithm for HYPERVOLUME in dimension d yields an algorithm for 
KMP in dimension \_d/2\ by Theorem 4, and since HYPERVOLUME is the easiest special case, we get the 
following corollary. 
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Corollary 6. // there is no improved algorithm for KMP, then there is no algorithm with runtime 
Q(^^yd/2\/2-e-^ /or any of Hypervolume, Unitcube-KMP, Cube-KMP, and fc-GROUNDED, for any e > 0. 

This shows a lower bound for all studied special cases. Note that there is, however, a wide gap to the best 
known upper bound of 0(n(''+^^/^) for HYPERVOLUME, Unitcube-KMP, and Cube-KMP. Also note that, 
other than the unconditional lower bound of fl{n\ogn) for HYPERVOLUME, there were no lower bounds for 
any of these special cases before. 

Furthermore, setting A; = 1, Theorem 4 immediately implies that the recent algorithm for 2-Grounded 
with runtime 0(n^''~^'/^ log^ n) [17] is optimal (apart from logarithmic factors and if there is no improved 
algorithm for KMP). 

Corollary 7. If there is no improved algorithm for KMP, then there is no algorithm for 2-Grounded with 

runtime 0{n^'^^^^/^^^) for any e > 0. 

An alternative way of stating this result would be that any algorithm faster than the one from [17] for 
2-Grounded would yield an improved algorithm for KMP. To simplify our runtime bounds, in some proofs 
we use the following technical lemma. Informally, it states that for any fc-GROUNDED algorithm with runtime 
T{n,d) we have T(0(n),d) ^ 0{T(n, d)). Intuitively, this should hold, since any reasonable algorithm for 
this problem has runtime T{n,d) = 0{n'^), i.e., bounded by a polynomial in n, and for polynomial T this 
simplification is valid. Note that in this paper we assume to be a constant. 

Lemma 8. Fix ^ k ^ d and c > 1. If there is an algorithm for ^-Grounded with runtime Tfc.GRouNDED('^, d) 
then there is another algorithm for A;-Grounded with runtime T'^.grounded('^> ^) satisfying 

Due to space constraints, the proofs of this and other lemmas can be found in the appendix. 
1.2 Notation and Organization 

A box is a set of the form B = [ai,6i] x ... x [0^,6^] C M'', ai,bi G R, ^ bi. A cube is a box with 
all side lengths equal, i.e., |6i — ai| = ... = \bd — ad\- Morever, a KMP instance is simply a set M of 
n boxes. In Cube-KMP all these boxes are cubes, and in Unitcube-KMP all these boxes are cubes of 
common side length. In HYPERVOLUME, all input boxes share the vertex (0, ...,0) e M**, i.e., each input 
box is of the form _B = [0, &i] x . . . x [0, bd]- In fc-GROUNDED, the projection of each input box to the first k 
dimensions is a Hypervolume box, meaning that each input box is of the form B = [ai,bi] x . . . x [od, bd] 
with ai = . . . = Ofc = 0. 

We write the usual Lebesgue measure of a set A C R"^ as vol{A). For sets i?, A C R'' we write vOLfj{A) := 
VOL(i? n A), the volume of A restricted to R. For a KMP instance M we let U{M) := [Jb^m ^- To shorten 
notation we write VOl(M) := V0l(ZY(M)) and VOLfl(M) := VOL(i? n i/(M)). 

In the next section we present the proof of Theorem 1. In Section 3 we prove Theorem 3. The proof of 
Theorem 4 is split into Section 4 and Section 5: We first give the reduction for 2- Grounded (again split into 
the case d = 1 and a generalization to larger dimensions) and then generalize this result to 2A;-Grounded, 
fc > 1. We close with an extensive list of open problems. 

2 Hypervolume ^ Unitcube-KMP 

In this section we prove Theorem 1 by giving a reduction from Hypervolume to Unitcube-KMP. This 
shows that the runtime for Hypervolume can be upper bounded by the runtime for Unitcube-KMP. 

We first describe the intuition, for clarification see Figure 1. Given an instance of Hypervolume, let A be 
the largest coordinate of any box. We extend all boxes to cubes of side length A, yielding a Unitcube-KMP 
instance. In this process, we make sure that the new parts of each box will not lie in the positive orthant 
M^Q, but in the other orthants, as depicted in Figure 1. This means that the restriction of each box to M^g 
stays the same, implying that the volume of the newly constructed cubes - restricted to R>q - is the same as 
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Figure 1: The left hand side depicts a Hypervolume instance consisting of 4 boxes in 2 dimensions. This instance 
is transformed into a Unitcube-KMP instance on the right hand side. We compute the volume of the constructed 
instance once with and once without the cube C, from which we can deduce the volume of the shaded part. The side 
length of C is A, the largest coordinate of any input box. 

the volume of the input boxes. To compute this restricted volume we introduce the cube C = [0, A]''. We 
compute the volume of the constructed Unitcube-KMP instance once with and once without the additional 
cube C. From this wc can infer the volume of the input HYPERVOLUME instance. In the following we give 
the details of this construction. 

Let M = {Bi, Bn\ be an instance of Hypervolume, i.e., Bi is a box in M*^ of the form 

Bi = [0,x[''>] X ... X [0,4'^]. 

Let A := ^- ' be the largest coordinate of the input instance. We construct two Unitcube-KMP 

instances M' and M^, where all cubes have side length A: The instance M'^ contains, for each 1 < i < n, the 
box 

B[:= [a.«-A,a=W]x...x[4)-A,4)]. 

We add the cube C := [0, A]'^ to get M' := U {C}. 

Note that the boxes C and are cubes of side length A and that we have CCiBl — Bi for each 1 ^ i ^ n. 
Hence, we have VOl(M) = volc(M^). We compute the latter using the inclusion-exclusion principle as 

voLc(AfO = vol(C) + voL(Af^) - vol(M'), 
which can be seen easily if we rephrase it giving more details, 

VOL(CnW(M^)) = vol(C) + vol(W(AO) - vol(CUW(MO) 
= vol(C) + vol(W(AO) - vol(W(M')). 

As VOl(C) = A'' is trivial, we reduced the input instance M to two Unitcube-KMP instances, M' and 
M^. Since these instances can be constructed in time 0{n), we showed that there is an algorithm for 
Hypervolume with runtime 

7hypervolume('^, d) < 2Tunitcube-KMp('^ + l,d) + 0{n). 

Since any algorithm for Unitcube-KMP at least has to read its complete input, implying 

T'uNiTouBE-KMp(n,c;) ^^{n), we get 

7hypervolume(^: d) ^ 0(7ijkitcube-KMP (n + l,d)). (1) 

We can slightly simplify this bound, decreasing n + 1 to n, by the technical Lemma 8. Plugging in the con- 
structed algorithm for Hypervolume with runtime THYPERvoLUME(f^j f^), this lemma guarantees an algorithm 
for Hypervolume with runtime 

^Hypervolume ('^j ^) ^ C'(2hypervolume(?^/2, d)) 
^ 0(T'hypervolume(?^ ~ ^)) 
(1) 

^ 0(Tunitoube- KMP {n, d)). 
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Figure 2: The left hand side depicts a Cube-KMP instance with smaUest cube C. The shaded region (with volume 
v) is the set of points contained in C but no other box. The result of our transformation is shown on the right hand 
side. All boxes now are cubes with C's side length, while inside C nothing changed. Computing the volume of the 
transformed instance once with and once without C yields v. 

3 Unitcube-KMP > Cube-KMP 

In this section we prove Theorem 3 by giving a reduction from CuBE-KMP to Unitcube-KMP. This shows 
that the runtimes of Cube-KMP and Unitcube-KMP arc polynomially related. 

To get an intuition, assume we are given any CuBE-KMP instance. We will show how to compute its 
volume assuming that we can solve Unitcube-KMP. For this, consider the smallest cube C in the input 
instanc;c, i.c;.. the one with smallest side length. We will compute the volume v of space that is contained 
in C but no other box of the input instance, see Figure 2. Then we recursively compute the volume w of the 
union of the n — 1 boxes that are left after deleting C. The sum of both, i; + u;, is the volume of the input 
instance. 

To compute v, we may modify each cube B ^ C in any manner such that the restriction of i? to C stays 
the same, this does not change v. We show that we can modify all other cubes such that their restriction 
to C remains untouched and they become cubes of C's side length. This is depicted in Figure 2. Applying 
this construction to all input boxes, we get a Unitcube-KMP instance that, inside C, looks the same as the 
input Cube-KMP instance. Now we compute the volume of this new instance, once with and once without 
the cube C. Subtracting both values, we get exactly v. This finishes the reduction. In the following we give 
the details of this construction. 

Let M = {Bi, . . . ,Bn} be an instance of CuBE-KMP. Sort the cubes Bi,. . . ,Bn descendingly by side 
length and consider the smallest cube C := B„. We want to compute the volume v of the space that is 
contained in C but no other box in M. With M* := M\ {C} we have 

V = vol(M) - vol(M*), 

since the extra volume in vol(M) compared to V0l(M*) stems from the space that is contained in box 
C — Bn but no other box Bj,. Once we know how to compute w, there is an immediate recursive algorithm 

to compute vol(M): 

1. Compute V. 

2. Recursively compute vol(M,). 

3. Return vol(M) = v + vol(M,). 

In the remainder of this section we describe how to reduce the computation of v to solving Unitcube- 
KMP instances. Write 

C = [aiM] X ... X [adM], 

and denote by A = foi — ai = . . . = 6^ — the side length of the cube C. Consider an input box Bi, i < n, 
that has non-empty intersection with C (all other boxes are unnecessary for computing v), with coordinates 

Bi = [xi,yi] X ... X [xd,yd]- 
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As Bi has larger or cqiial side length compared to C, we cannot have < Xi ^ yi < hi for any i. Instead, 
for each dimension i, there are the following three possible relations between the intervals [oj, bi] and [xi,yi]. 
In each case we define an interval S- of length A satisfying B'^ n [oj, bi] = [xi,yi] fl [ui, bi]: 

[xi,Xi + A]. 

[y^ - A,yi]. 

[ai,bi]. 

Then setting B' = B[ x . . . x B'^ wc get a cube of side length A with B' (IC = BdC. This completes the 
construction of a modified cube B' that is of C's side length and looks the same as B when restricted to C. 

Now we construct two Unitcube-KMP instances M' and as follows. M' contains the modified 
version of every input cube -Bj, 1 < i ^ n, that has non-empty intersection with C. In particular, M' 
contains C, since C = C. Moreover, we set = M' \ {C}. Now consider 

vol(M')-vol(M^). 

This measures the volume of space contained in C but no box in . As the restriction of every box in 
to C is the same as the restriction of the corresponding box in M* to C, we get 

vol(M') - vol(M^) = vol(M) - vol(M*) = v. 

Thus, we can indeed compute v by solving two Unitcube-KMP instances, which completes the reduction. 

What is the runtime of the constructed algorithm? Using the above recursive procedure, we reduce 
the computation of Cube-KMP to 0{n) instances of Unitcube-KMP. After we initially sort the boxes, 
constructing one of these instances takes time 0{n). Thus, in total we need additional time 0{n^), yielding 
a runtime bound of 

7cuBE-KMp(n, d) < 0{n ■ TuNiTcuBE-KMp(n, d)) + 0(n^). 

As any algorithm for Unitcube-KMP at least needs to read its whole input, we have Iunitcube-kmp d) ^ 
Cl{n), which allows to further simplify the above time bound to the desired bound 

TcvBE-KMp{n, d) ^ 0(n ■ Tunitcube-kmp {n,d)). 

4 2-Grounded ^ KMP 

We first show the reduction of Theorem 4 for 2-Grounded, i.e., we show TuMpin, d) < 0{T2-GRo\jNDED{n,d + 
1)) by giving a reduction from KMP to 2-Grounded. This already implies Corollary 7 and lays the foun- 
dations for the complete reduction given in the next section. 

We begin by showing the reduction for d = 1. As a second step we show how to generalize this to larger 
dimensions. 

4.1 Dimension d = 1 

We want to give a reduction from KMP in 1 dimension to 2-Grounded in 2 dimensions. Note that the 
latter is the same as Hypervolume in 2 dimensions. Let M be an instance of KMP in 1 dimension, i.e., a 
set of n intervals in R. We will reduce the computation of vol(M) to two instances of 2-Grounded. 

Denote by xi < . . . < Xm the endpoints of all intervals in M (if all endpoints are distinct then m = 2n). 
We can assume that a;i = after translation. Consider the boxes 

Ai := [m - i - l,m - i] X [xj, Xj+i] 

in for 1 ^ ? ^ m — 1, as depicted in Figure 3. Denote the union of these boxes by A. Note that the 
volume of box Ai is the same as the length of the interval [xj,Xj-)-i]. This means that we took the chain of 



• If fli < Xi ^ 5i ^ j/ii Set B[ = 

• li Xi ^ tti ^ yi < bi'. Set = 

• If Xj < ttj < 6i < yii Set B- = 



6 



XI X2 X3 X4X5 



X6 



X5 

xn 



As 








A4 




Ci 




A3 










A2 




1 


1 1 1 


— 



Figure 3: The left hand side depicts all endpoints = xi < . . . < a;6 of a 1-dimensional KMP instance. An input 
interval I is indicated. The right hand side shows the result of our transformation. Each interval [xj, Si+i] to the left 
corresponds to a box Ai to the right. The interval I gets mapped to the box C/. The shaded regions depict the set 
A (■, the union of all Ai) and the set To ( ). 



intervals {[xi, Xj+i]} and made it into a staircase of boxes {Ai}, where each box has the same volume as the 
corresponding interval. 

Now consider an interval / = [xj,Xk] € M. We construct the box 

Ci := [0,m-j] X [0,a;fe], 

also shown in Figure 3. Then C/ contains the boxes Ai with j ^ i < k and (its interior) has no common 
intersection with any other box Ai. This is easily seen asAiCCiiSm — i^m~j and x^+i ^ Xk- Hence, 
for any interval / G M we constructed a box C/ that contains exactly those boxes Ai whose corresponding 
interval [a;,,a;i+i] is contained in /, or in other words 

[xi, Xi+i] C / <S4> C C/, 
voL([a;i, a;i+i] n 7) = voL{Ai n d). 

From these properties it follows that the volume of C/ restricted to A is the same as the length of I, i.e., 

vola(C/) = vol(/). 

Furthermore, considering the whole set M of intervals, the interval is contained in some interval 

in M iff the box Ai is contained in some box in Cm '■= {Cj \ I G M}. This yields 

VOL(M) = VOLa(Cm)- 

It remains to reduce the computation of VOLa (Cm) to two 2-Grounded instances. For this we consider 

To := [J C[^i,xi]- 

Informally speaking. To consists of all points "below" A, as depicted in Figure 3. Note that no set Aj is 
contained in Tq. Moreover, we consider the set Ti := Tq U A. Observe that we can write 

Ti = IJ C[^i,xi+i], 

l^i^rn—l 

since Ai C C[xi,xi+i]- Note that both sets Tq and Ti are unions of 0{n) 2-Grounded boxes. Informally, Tq is 
the maximum 2-Grounded instance that has VOL^(To) = 0, and Ti is the minimum 2-Grounded instance 
with vola(Ti) = vol(^). Now, we can compute VOLa(Cm) as follows. 

Lemma 9. In the above situation we have 

vol^(Cm) = vol(A) +vol(To UW(Cm)) - vol(Ti UW(Cm)). 



7 



Proof. Set U := U{Cm)- Using Tq C Ti and ^ = Ti \ Tq in a sequence of simple transformations, we get 

vol(Ti U [/) - vol(To U [/) = vol((Ti V^U)\ (To U U)) 

= VOL((Ti \ To) \ U) 

= \oh{A \ U) 

= vol(^) - vol{A n U) 

= vol(A) -vol^ (Cm), 

which proves the claim. □ 

Note that vol(A) = X^jVOL(Ai) — J^i ^ Xi\ = \xm — xi\ is trivial. Also note that both sets Tg 
and Ti are the union of 0{n) 2-Grounded boxes, so that VOL(rb \JU{Cm)) can be seen as a 2-Grounded 
instance of size 0{n), for both b G {0, 1}. Hence, we reduced the computation of the input instance's volume 
VOl(M) to V0L/i(Cm) and further to the 2-Grounded instances VOL(To \JU{Cm)) and VOL(Ti \JU{Cm))- 

As we have to sort the given intervals first, we get 

rKMp(n, 1) ^ C(r2.GRouNDED(0(n),2) + nlogn). 

Note that this inequality alone gives no new information, as already Klee [13] showed that Tkmp(?^, 1) ^ 
©(nlogn). However, we get interesting results when we generalize this reduction to higher dimensions in the 
next section. 



4.2 Larger Dimensions 

In this section wc show that the reduction from the last section easily carries over to larger dimensions, yielding 
a reduction from KMP in d dimensions to 2-Grounded in d + 1 dimensions. This implies TKMv{n,d) < 
C(^2-Grounded(?^, ci+ !))• However, here we do not give all details of the proof, since we prove a more general 
statement in the next section. 

Assume we are given a KMP instance M in dimension d. The idea is that we use the dimension doubling 
reduction from the last section on the first dimension and leave all other dimensions untouched. More 
precisely, for a box B d M let i^i{B) be its projection onto the first dimension and let 7r*(i3) be its projection 
onto the last d — 1 dimensions, so that B = 7ri{B) x Tr^,{B). Now follow the reduction from the last section 
on the instance M' := {tti{B) \ B e M}. This yields sets A, Tq, Ti, and a box C/ for each / G M'. 

We set Cb '■= C^^^i^b) x 7r*(-B) and Cm = {Cb \ B e M}. A possible way of generalizing A would be to set 
A" := Ax Then we would be interested in VOhA"{CM), which can be seen to be exactly vol(M). This 

definition of A" is, however, not simple enough, as it is not a difference of 2-Grounded instances (unlike 
A = Ti\To). To give a different definition, assume (after translation) that all coordinates of the input instance 
are nonnegative and let A be the maximal coordinate in any dimension. We set A' := Ax [0, A]''^^ and still 
get the same volume VOLa'(Cm) = VOl(M). This allows to generalize Tq and Ti to Tq := Tq x [0, A]''"^ and 
T{ := Ti X [0, A]'^-S while still having 

voLA'ic&i) = vol(A') + voL(r^ uu{Cm)) - VOL(r{ U U(Cm)). 

Note that Tq and T{ are also a union of 0{n) 2-Grounded boxes, so a volume such as vol(Tq U U{Cm)) 
can be seen as a 2-Grounded instance. This completes the reduction and yields the time bound 

T'KMp(n,d) 0(r2.GRouNDED(C(n),(i + 1) + nlogn). 

Using the lower bound 12 (nlogn) of Corollary 2 we can hide the additional nlogn in the first summand. 
We may simplify further by the technical Lemma 8: By using the algorithm guaranteed by Lemma 8 with 
runtime T^.q^q^ndedC'^' ^) S^t a runtime bound of 

TKMp(n,d) SS 0(r^.GEOUNDED(C?W,C?+ 1)) 0(T2.GR0UNDED(n,d+ 1)), 

which implies Corollary 7. 
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5 2k-Grounded ^ KMP 

In this section we now prove the full version of Theorem 4, i.e., we give a reduction from KMP in dimension d 
to 2fc-GROUNDED in dimension d+k. For this we enrich the reduction from the last section with an advanced 
application of the inclusion-exclusion principle. 

Let M be a KMP instance. After translation we can assume that in every dimension the minimal 
coordinate among all boxes in M is 0. Denoting the largest coordinate of any box by A we thus have 
B C [0, A]'' = n'^ for all B € M, where n := [0, A]. 

The first steps of generalizing the reduction from the last section to the case fc > 1 are straight forward. 
We want to use the dimension doubling reduction from Section 4.1 on each one of the first k dimensions. For 
any box B £ R''- denote its projection onto the i-th dimension by ni{B), 1 ^ i ^ k, and its projection onto 
dimensions fc + 1, . . . , d by 7r*(i3). We use the reduction from Section 4.1 on each dimension 1 < i ^ fc, i.e., 
on each instance := {77^(5) | B G M}, yielding sets A^'\T^'\t^\ and a box c|'^ for each / e M^. 
We assume that all these sets are contained in [0, A]^ = fl^, meaning that all coordinates are upper bounded 
by A (this holds after possibly increasing the A we chose before). 

For a box B £ M we now define 

GB:=Ci%x...xCi';^^^xMB). 

This is a box in K'^+'=, it is even a 2A;-Grounded box, as its projection onto the first 2k coordinates has the 
vertex (0, . . . ,0). Set Cm '■= {Cb \ B e M}. As shown by the following lemma, we want to determine the 
volume VOLa (Cm) in 

A:= A(i) X ... x^^'') X [0,A]'^-'=. 
Lemma 10. In the above situation we have 

vol(M) = vola (Cm)- 

Unfortunately, the set A is not simply a difference of two 2fc-GROUNDED instances. Thus, the hard part 
is to reduce the computation of VOLa(Cm) to 2A;-Grounded instances, which we will do in the remainder 
of this section. For 1 < i < fc and b G {0, 1} we set 

f^W := X T^^ X n-^+^-^K 

This set in $1'^+*^ consists of all points x whose projection to dimensions 2i — 1 and 2i is contained in T^J'"' . Note 
that each set T^^'^ can be written as the union of 0(n) 2A;-Grounded boxes, since T^^^ is the union of ^ = 0{n) 
2-Grounded boxes in R?, i.e., T^'' = {f-^^ Cj, so that we may write f^^ = (J^^i ^^('-i) x Cj x 

Thus, we can use an algorithm for 2fc-GROUNDED to compute any volume of the form YOl.{T^^^ \JV) , where V 
is a union of 0{n) 2fc-GROUNDED boxes. 
Furthermore, define for S C [A;] 

^ies ^ ie[k]\s 

Note that Ds Q Ds' holds for S C S'. We can express A using the sets Ds as shown by the following lemma. 
Lemma 11. In the above situation we have 

A= f] D{^\D^. 

Moreover, each Ds can be written as the union of 0{n) 2A;-Grounded instances, since the same was true 
for the sets Tj|*^ Hence, we can use an algorithm for 2fc-GROUNDED to compute the volume 

Hs := VOL(£)sUW(Cm)). 

Next we show that we can compute VOL^ (Cm) from the Hs by an interesting usage of the inclusion-exclusion 
principle. 
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Lemma 12. In the above situation we have 

vola(Cm) = VOL(A) + J2 (-l)'^'^^s. 

SC[fe] 

As vol(A) = A''"'^ • Hi^j^fc vol(j4(')) is trivial, we have reduced the computation of VOLa (Cm) to 

2k _ 0(1) instances of 2fc-GROUNDED, each consisting of 0{n) boxes. During the construction of these 
instances we need to sort the coordinates, so that we need additional time 0{nlogn). This yields 

TKMp{n, d) ^ C(r2fe-GRouNDED(C(n), d + k) + nlogn). 

Because of the lower bound from Corollary 2, we have T2k-GRovNDEB{0{n),d + k) = O(nlogn), so we can hide 
the second summand in the first, 

TKMp{n,d) < 0{T2k -Grounded {0{n),d + k)). 

We may use the technical Lemma 8 to get rid of the inner O: This lemma guarantees an algorithm with 
runtime T^k-GRovNOEoii^^d) such that we get 

-Grounded 

(0(n), d + k))^ 0(T2fe.GRouNDED(n, d + k)). 

This finishes the proof. 

6 Conclusion 

We presented reductions between the special cases Cube-KMP, Unitcube-KMP, Hypervolume, and k- 
Grounded of Klee's measure problem. These reductions imply statements about the runtime needed for 
these problem variants. We established Hypervolume as the easiest among all studied special cases, and 
showed that the variants Cube-KMP and Unitcube-KMP have polynomially related runtimes. Moreover, 
we presented a reduction from the general case of KMP to 2A;-Grounded. This allows to transfer W[l]- 
hardness from KMP to all special cases, proving that no n°('') algorithm exists for any of the special cases 
under reasonable complexity theoretic assumptions. Moreover, assuming that no improved algorithm exists 
for KMP, we get a tight lower bound for a recent algorithm for 2-Grounded, and a lower bound of roughly 
j^(<i-i)/4 fQj. a,ii other special cases. Thus, we established some order among the special cases of Klee's measure 
problem. 

Our results lead to a number of open problems, both asking for new upper and lower bounds: 

• Is there a polynomial relation between Hypervolume and Unitcube-KMP, similar to the relation 
between Cube-KMP and Unitcube-KMP, or do both problems have significantly different runtimes? 

• Show that no improved algorithm exists for KMP, e.g., assuming the Strong Exponential Time Hy- 
pothesis, as has been done for the Dominating Set problem, see [14]. Or give an improved algorithm. 

• Assuming that no improved algorithm for KMP exists, we know that the optimal runtimes of Hyper- 
volume and Cube-KMP/Unitcube-KMP are of the form n'^-d+o{d)^ .^j^j^ ^ g ^^^^^ ]^^3j Determine 

the correct value of c. 

• Generalize the ©(n^''"^)/^ log^ n) algorithm for 2-Grounded [17] to an 0{n'-'^~''^^^~^°'-^^ ) algorithm for 
2fc-GROUNDED. This would again be optimal by Theorem 4. 

• We showed the relation TKMp{n,d) ^ O(T2fe_GR0UNDED('T-, + k)). Show an inequality in the opposite 
direction, i.e., a statement of the form Tfc.GRouNDED('^, c^) ^ 0{TKMp{n,d')) with d' < d. 
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A Proof of Lemma 8 



Lemma 8. Fix ^ k ^ d and c > 1. If there is an algorithm for fc-GROUNDED with runtime T/j.qroundedI^^, d) 
then there is another algorithm for A;-Grounded with runtime T'^.qroundedC'^' satisfying 

Proof. Let M be an instance of A;-Grounded of size \M\ = n. Denote by < . . . < Z2n the coordinates 
in the first dimension of all boxes in M. Let a\ := Z(i_a)n and bi := where a := l/(3d). Denote 

by {xi < ai} the set of all points x = {xi, . . . ,Xci) G M'' with xi < ai, similarly for {ai < a;i < 61} and 
{xi > 61}. We consider the three KMP instances 

Ma := {B n {xi ai} I B e M}, 
Mb := {B n {xi ^bi}\B & M}, 
M' := {B n {ai ^ Xi ^ bi} \ B e M}. 

Note that all three instances can be seen as fc-GROUNDED instances: Projected onto the first k dimensions, 
all boxes in Ma share the vertex (0, . . . , 0), all boxes in M{, share the vertex (61 , 0, . . . , 0), and all boxes in M' 
share the vertex (ai, 0, . . . , 0), so after translation they all share the vertex (0, . . . , 0). Moreover, {xi < 04} 
and {xi > 61} contain only (1 — a)n coordinates of boxes in M in the first dimension. Hence, there are at 
most (1 — a)n boxes intersecting {xi < ai}, so after deleting boxes with volume we get \Ma\ ^ (1 — ce)n, 
similarly for Mf,. This reasoning docs not work for M', it might even be that all n boxes arc; present in M': 
If a box has left coordinate smaller than ai and right coordinate larger than bi , then none of its coordinates 
is seen in {ai ^ 0:1 ^ fei}, although it has non-empty intersection with {oi ^ xi ^ 61}. However, such a box 
in M' is trivial in the first dinicmsion: its coordinates in the first dimension are simply [ai, 61]. If all boxes in 
M' were trivial in the first dimension, then M' would clearly be simpler than the input instance. Although 
this is not the case, we can bound the number of boxes in M' that are non-trivial in the first dimension: 
Since there are at most 2an coordinates z,; in [oi, 61], all but at most 2an boxes in M' are trivial in the first 
dimension. Thus, also M' is easier than the input instance M , in a certain sense. 

Note that VOl(M) = vOL(Ma) + VOL(Mf,) + vol(M') and Ma, Mb are strictly easier than M, as they 
contain at most (1 — a)n boxes. We have to simplify M' further. For this, we use the same construction as 
above (on M and dimension 1) on M' and dimension 2, i.e., we split by coordinates in dimension 2 at 02 
and 62. This yields three fc-GROUNDED instances. Two of them contain at most (1 — a)n boxes. The third 
one, M", may contain up to n boxes. However, all but at most 4an of these boxes are trivial in the first and 
second dimension, meaning that their projection onto the first 2 dimensions is [ai,6i] x [02,62]- 

Iterating this reduction d times yields 2d instances of fc-GROUNDED containing at most (1 — a)n points 
and one instance M* that may contain up to n boxes. However, all but at most 2d an of these boxes are 
trivial in all d dimensions, meaning that they are equal to [oi, 61] x ... x [a^, bd\. Since all boxes in M* are 
contained in [ai, &i] x . . . x [aj,, bj], if any such trivial box exists, the volume of M* is trivial. Otherwise M* 
only contains at most 2dan = |n < (1 — l/{3d))n = (1 — a)n boxes. Thus, we have reduced the computation 
of VOL(M) to 2d + 1 instances of fc-GROUNDED with at most (1 — a)n boxes each. The reduction itself 
can be made to run in 0{n) time. Hence, if we solve the reduced problems by an algorithm with runtime 
T/c-GRouNDED(n, d), then we get an algorithm with runtime TX-groundedI'^' ^) satisfying 

^fc-GROUNDED d) ^ (2rf + l)Ifc.GROUNDED((l " d) + 0{n). 

As every algorithm for fc-GROUNDED at least has to read its whole input, we can hide the 0{n) by the first 
term, 

^fc-GROUNDED ('^i'^) < C(Ife-GROUNDED((l " a)n,d)), 

or, 

^fe-GROUNDEDW(l -a),d) ^ 0(Tfe.GROUNDED(n,d)). 
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Repeating this construction an appropriate number of times we can increase the constant 1/(1 — a) to any 
constant c > 1, while the factor on the right hand side is stih bounded by a constant. This finally yields an 
algorithm with runtime satisfying 

-^fc-GROUNDED ^ C(Tfe-GROUNDED('T') f^))- 

□ 

B Proof of Lemma 10 

Lemma 10. In the situation of Section 5 we have 

vol(M) = vola(Cm)- 

Proof. Denote by a;^*^ < . . . < the coordinates of all boxes in M in the i-th dimension. We can express 
VOL(M) in terms of the boxes 

Eji:-,3d '■= [x\^iX^jXi^ X ... X ,a;y^_J, 

for 1 < ji < TOj. Since each such box is either completely included in some box in M or does not contribute 
to vol(M), we have 

VOL(M)= ^ CiY(M)]-voL(£;,„...,,J, 

where [X] is 1 if X is true, and otherwise. 

Recall from the reduction in Section 4.1 that ^^'^ = A^'^ U. . .UA^^._^, and there is a one-to-one correspon- 
dence between intervals [x^*\a;j!^]^] and 2-dimcnsional boxes A^^\ in particular both have the same volume. 
This carries over to a one-to-one correspondence between d-dimensional boxes Eji,...,j^ and (d-l-fc)-dimensional 
boxes 

in particular both have the same volume. 

Additionally, recall that an interval / e M*^'^ includes [a;^*\ x^^^j^] if and only if the 2-dimcnsional box cj*'' 

contains A^'''. If / does not include [x'j\x'jl-^^\, then vol(/ n [x'j'\ x'^^^l^]) = 0, and we also have vOL(cj'^ n 

Aj^^) = 0. Hence, we have for any B € M that Ej^^,,,j^ C B if and only if E'^^ ^^ C Cb, which implies that 
C U{M) if and only if Ej^^ j^ C U{Cm)- Furthermore, E'j^ j^ is either included in ^/(Cm) or does 
not contribute to vol(Cm)- In total, we get 

VOL(M)= [En,-,jd^^{M)]-VOL{Ej,,...j,) 

ji,---,jd 

= E [£^;,...j.cZ^(Cm)]-vol(£;;,...,,J = vol^(Cm), 
ji,---,jd 

which completes the proof. □ 



C Proof of Lemma 11 

Lemma 11. In the situation of Section 5 we have 

A= f] D[iy\D<i^. 
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Proof. We have = Ui£[fe] Tq'^ and f^'^ C f^''^ for all i, implying 



This yields 



l^i^fe ^l^i^fc ^ ie[k] 

A point a; = (a;i, . . . , Xd+fe) e is in the set on the right hand side if and only if it has the following three 

properties: 

• Xi € [0, A] for all 1 ^ i < d + k, 

• {x2i-i, X2i) lies in T^^ for all 1 < i < A;, 

• {x2i-i,X2i) does not lie in Tq^ for all 1 < i < fc. 

Since A^'^ = T^^ \ T^\ this description captures exactly A^^) x ... x A^''^ x [0,A]'^-'= = A, finishing the 
proof. □ 

D Proof of Lemma 12 

Lemma 12. In the situation of Section 5 we have 

vola(Cm) = VOL(A) + Y (-l)'^'^^s. 

SC[fe] 

Proof. In this proof we write for short U :=U{Cm)- We clearly have 

VOL(^) - VOLA(f/) = V0L(A \ U). (2) 

We first show 

VOL(^\l7)= Y (-l)'^' + '(^5-i^0), (3) 

07^SC[fc] 

and simplify the right hand side later. Using A = Hi^i^fe ^{i} \ ^0 (Lemma 11) and the inclusion-exclusion 
principle we arrive at 

VOL(A\i7) = VOL( fl D^ijMD^UU)) 

= Y {-lf^^'yOL([jD{,y\{D,UU)). 

Note that Ds = Uies ^{i} ) so that the above equation simplifies to 



VOL(A\[/)= Y (-1)""+Vol(i)s\(D0U;7)) (4) 

0#SC[fc] 



Using the definition of Hg and C for any S C [k], we get 

Hs- HiD= vol{Ds UU) - voL(i:)0 U U) 
= VOLiiDs U U) \ [D^ U [/)) 
= yoL{Ds \ {D^ U U)). 
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Plugging this into (4) yields (3). 
Observe that 

05^SC[fe] 

This allows to further simplify (3) to 

voL{A\U)= 5^ (-1)1 

SC[fe] 

Plugging this into (2) yields the desired equation. 
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